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Abstract 

Precision comparisons of different atomic frequency standards over a period of a 
few years can be used for a sensitive search for temporal variations of fundamental 
constants. We present recent frequency measurements of the 688 THz transition in 
the ^^^Yb+ ion. For this transition frequency a record over six years is now available, 
showing that a possible frequency drift relative to a cesium clock can be constrained to 
(—0.54 ± 0.97) Hz/yr, i.e. at the level of 2 • 10~^^ per year. Combined with precision 
frequency measurements of an optical frequency in ^^^Hg"^ and of the hyperfine ground 
state splitting in *''Rb a stringent limit on temporal variations of the fine structure 
constant a: dlna/dt — (—0.26 ± 0.39) • 10~^^ yr~^ and a model-dependent limit for 
variations of the proton-to-electron mass ratio fi in the present epoch can be derived: 
dlnfj,/dt = ( — 1.2 ± 2.2) ■ 10"^'' yr~^. We discuss these results in the context of astro- 
physical observations that apparently indicate changes in both of these constants over 
the last 5-10 billion years. 



1. Search for variations of constants in astrophysics and with atomic clocks 

Over the past few years there has been great interest in the possibihty that the funda- 
mental constants of nature might show temporal variations over cosmological time scales 
[U m [31 m [S] . Such an effect - as incompatible as it seems with the present foundations of 
physics - appears quite naturally in the attempt to find a unified theory of the fundamental 
interactions. An active search for an indication of variable constants is pursued mainly in 
two areas: observational astrophysics and laboratory experiments with atomic frequency 
standards. The two most important quantities under test are Sommerfeld's fine structure 
constant a — e'^ /{Aneotic) and the proton-to-clcctron mass ratio /i — nip/me. While a is the 
coupling constant of the electromagnetic interaction, fi also depends on the strength of the 
strong force and on the quark masses via the proton mass. It is important that these quan- 
tities are dimensionless numbers so that results can be interpreted independently from the 
conventions of a specific system of units. Quite conveniently, both constants appear promi- 
nently in atomic and molecular transition energies: a in atomic fine structure splittings and 
other relativistic contributions, and /i in molecular vibration and rotation frequencies as well 
as in hyperfine structure. 

A multitude of data on variations of a and fi has been obtained from astrophysical 
observations but the present picture that is obtained is not completely consistent: Evidence 
for a variation of a has been derived from a shift of wavelengths of metal ion absorption 
lines produced by interstellar clouds in the light from distant quasars [B] . These observations 
suggest that about 10 billion years ago (redshift range 0.5 < z < 3.5), the value of a was 
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smaller than today by Aa/a = (—0.543 ± 0.116) • 10~^, representing A.7a evidence for a 
varying a [7]. (Here the definition is used that a negative sign of Aa corresponds to a 
temporal increase of the constant.) Assuming a linear increase of a with time, this would 
correspond to a drift rate d\na/dt = (6.40 ± 1.35) • 10~^^ yr~^ [7]. Other evaluations of 
spectra obtained with a different telescope and using other selections of quasar absorption 
systems reach similar or even higher sensitivity for Aa/a but are consistent with Aa = 
for all look-back times [H [9] . 

In the case of the proton-to-electron mass ratio a relative change A/i//i = (2.4±0.6) -10"^ 
has been inferred from the analysis of absorption systems of H2 clouds at redshifts z = 2.6 
and 3.0 (look-back time w 12 Gyr), observed in the spectra of two quasars [TD]. Assuming a 
linear time evolution, this would correspond to dln^/dt « —2 • 10~^^ yr""'^. Looking only at 
the positive reports about variations, one would conclude \ Afj,/fj,\ > |Aa/Q;|. Indeed, in the 
framework of theories of unification it is plausible that the scale of the strong interaction 
shows a bigger rate of change than the electromagnetic force [TT| \T2\ [T3] . 

Other recent reports derive stringent limits on variations: an analysis of 18 cm OH lines 
at red shift z = 0.765 (look-back time 6.7 billion years) constrains the variation of the 
quantity F = ^^(a^^)!-" (g^: g-factor of the proton) to AF/F = (0.44 ± 1.07) • 10"^ [Ti] . 
A combination of the H ground state hyperfinc splitting interval (21 cm line) and metal 
ion absorption lines in quasar spectra can be used to test the constancy of the parameter 
X = a'^gp/fj, and a result AX/X = (1.17 ± 1.01) ■ 10~^ has been obtained for the redshift 
range 0.2 < z < 2.0. [15J Combined with either the non-zero "7] or the zero [3] value for Aa/a 
one obtains A/i//i = (2.31 ± 1.03) • 10"^ or A^l/^I = (1.29 ± 1.01) • 10"^, respectively. It 
should be kept in mind that the uncertainties cited here are usually the result of averaging 
over a large ensemble of data points with much higher individual statistical uncertainties. 
The significance of the results that detect non-zero variations therefore hinges on a complete 
understanding of all systematic effects that may introduce bias or correlations in the data 
and that may not be easy to detect. 

The approach of using precision laboratory experiments with frequency standards to 
search for variations of the fundamental constants [16| has the advantage that the relevant 
parameters are under the control of the experimenter, permitting a detailed investigation 
of possible systematic effects. The obvious disadvantage is that the limited duration of the 
experiments of typically a few years only is not well adapted to search for evolution on a 
cosmological time scale. Application of the methods of laser cooling and trapping has led to 
significant improvements in the precision of atomic clocks and frequency standards over the 
last years [ITl [18] : Primary cesium clocks based on atomic fountains [HI [20] realize the SI 
second with a relative uncertainty below 1 • 10~^^ and a similar level of accuracy has been 
achieved for the ground state hyperfine frequency of *^Rb [20]. Research towards optical 
clocks with trapped ions and atoms [HJ [TS] has resulted in several systems that approach 
the accuracy of primary cesium clocks or show an even higher reproducibility. In typical 
experiments performed so far to search for temporal changes of fundamental constants, 
an atomic transition frequency /at is measured with respect to a cesium clock (i.e. the 
frequency ratio /at//cs is determined, where fcs = 9192 631770 Hz, the SI value of the 
ground state hyperfine splitting frequency of ^^'^Cs). One measurement usually takes a few 
days of averaging time in order to reach a low statistical uncertainty. After a few years the 
measurement is repeated and changes in the frequency ratio can be looked for. 

This procedure has recently been undertaken for: the ground state hyperfine frequency 
of ®^Rb [22l[20], the 15 ^ 25 two-photon transition in atomic hydrogen [23l, the transition 
^51/2 ^i?5/2 at 1065 THz in Hg+ O [25], and the ^51/2 ^ ^£>3/2 transition at 688 THz 
in ^^^Yb+ (this work). Significant changes of a or /i would produce a clear signature in 
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these experiments, because they would influence the transition frequencies differently [Hl26j. 
Atomic gross structure scales with the Rydberg constant i?oo, hyperfine structure with the 
product of a^Roo and the nuclear magnetic moment. In addition, the fine structure constant 
a appears in relativistic contributions to the level energies that increase with the square of 
the nuclear charge and are consequently much more important for heavy atoms '27\. The 
sensitivity of a specific transition frequency to changes of a can be derived from ah initio 
relativistic atomic structure calculations [28l [29] . While measurements on a single atomic 
gross structure transition relative to a cesium clock can test variations of a combination 
of a and the cesium magnetic moment only, a combination of data from different optical 
transition frequencies can be used to obtain a model-independent result about variations 
of a alone [H [531 [SOj ■ The published most stringent laboratory limit on variations of the 
fine structure constant in the present epoch had been obtained by this method in 2004: 
dhia/dt — (—0.3 ± 2.0) • 10"^^ yr~^ [30 , using the three optical transitions mentioned 
above. 

Contrary to the astrophysical observations, all present laboratory experiments are com- 
patible with constancy of constants to within Icr. The question arises to which degree the 
results of the quasar absorption spectra and the laboratory experiments can be interrelated, 
because they investigate very different time scales and different regions of the universe. It 
has been pointed out recently that cosmological variations of physical constants linked to 
scalar fields should also be detectable on the surface of local gravitationally bound systems 
[31j . This means that terrestrial laboratory experiments should in principle be able to reflect 
variations observed in interstellar gas clouds. Concerning the possible time dependence, one 
may expect the most dramatic variations during the very early phases of the evolution of 
the universe, e.g. during the inflationary period. Because the time scale on which variations 
seem to appear in the astrophysical observations is logarithmically much closer to our epoch 
than to the Big Bang, a simple linear time evolution of the constants may be used as a first 
approximation to model the behavior over the past few billion years. 

In the following, we present recent precision measurements with the Yb+ optical fre- 
quency standard at PTB that augment the existing record of data from this system. In 
combination with data from NIST, Boulder, on the Hg^ standard a new model-independent 
limit on the variation of a is obtained. Combining this with data from LNE-SYRTE, Paris, 
on the constancy of the Rb/Cs hyperfine frequency ratio and using the Schmidt model for 
the nuclear g-factors of Rb and Cs, a limit on the variation of n is obtained. 

2. Measurements of the Yb+ optical frequency standard at 688 THz 

The ^^^Yb"'' ion is attractive as an optical frequency standard because it offers narrow 
reference transitions with small systematic frequency shift. At PTB, the electric quadrupole 
transition {^Sii2,F = 0) -> (^Dg/a, F = 2) at 436 nm wavelength (688 THz frequency) with 
a natural linewidth of 3.1 Hz is investigated [32l[33l[34j[35]. A single ytterbium ion is trapped 
in a miniature Paul trap and is laser-cooled to a sub-millikelvin temperature by exciting the 
low-frequency wing of the quasi-cyclic {F = 1) ^ {F = 0) component of the — > 
resonance transition at 370 nm (cf. Fig. 1). The S D electric quadrupole reference 
transition is probed by the frequency doubled radiation from a diode laser emitting at 871 
nm. The short term frequency stability of this laser is derived from a temperature-stabilized 
and seismically isolated high-finesse reference cavity. 

Figure lb shows a high-resolution excitation spectrum obtained with 90 ms long laser 
pulses, leading to an approximately Fourier-limited linewidth of 10 Hz, or a resolution A1//1/ 
of 1.4 • 10~^^. Since the duration of the probe pulse is longer than the lifetime of the excited 
state (51 ms), the observed maximum excitation probability is limited by spontaneous decay. 
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Figure 1: a) Simplified level scheme of the ^^^Yb+ ion. b) Excitation spectrum of the 
'^Si/2{F — 0,mF = 0) ^ ^Dt,/2{F = 2,mF = 0) transition, obtained with 20 probe cycles 
for each value of the detuning. The lincwidth of 10 Hz is approximately at the Fourier limit 
for the employed probe pulses of 90 ms duration. 



In order to operate the system as a frequency standard, both wings of the resonance are 
probed alternately, and the probe light frequency is stabilized to the line center according 
to the difference of the measured excitation probabilities ^34j . 

For quantitative studies of systematic frequency shifts we have compared the line center 
frequencies of two ^^'^Yb"'" ions stored in separate traps [33 . The mean frequency difference 
between the two traps was in the sub-hertz range for nominally unperturbed conditions. 
The dominant source of systematic uncertainty in the ^^^Yb+ optical frequency standard is 
the so-called quadrupole shift of the atomic transition frequency. It is due to the interaction 
of the electric quadrupole moment of the -D3/2 state with the gradient of static electric patch 
fields in the trap. The frequency difference between the two systems shows a relative insta- 
bility (Allan deviation) of crj,(1000s) = 5 • 10~^^, limited essentially by quantum projection 
noise [M] . 

The absolute frequency of the ^"^^Yb^ standard at 688 THz was measured relative to 
PTE's cesium fountain clock CSFl [SGlIS^. The link between optical and microwave fre- 
quencies is established with a femtosecond-laser frequency comb generator [38]. The ^^^Yb+ 
frequency was measured first in December 2000 (3^. Fig. 2 shows the results of all frequency 
measurements performed up to June 2006, showing a steady decrease of the measurement 
uncertainty and excellent overall consistency. Technical improvements over the years were 
made in the resolution of the ionic resonance, in the control of trapping conditions, and in 
the frequency comb generator setup. The latter was based on a titanium sapphire laser in 
the earlier experiments (until MJD 53000) and on an erbium doped fiber laser since [3D]. 
Between July 2005 and June 2006 five absolute frequency measurements with continuous 
averaging times of up to 36 h were performed. The obtained statistical uncertainties (about 
6 • 10^^^ after two days) are dominated by the white frequency noise of the cesium fountain. 
The relative systematic uncertainty contribution of the cesium reference of 1.2-10^^^ is based 
on a recently completed evaluation of CSFl that is applicable to the 2005/06 measurements 
in retrospect. The systematic uncertainty contribution for the Yb+ standard of 1.5 • 10""'^^ 
is largely dominated by an estimate of 1 Hz for the stray-field induced quadrupole shift. 

The weighted average over these measurements gives the present result for the frequency 
of the "iYb+ 25i/2(F = 0) ^ '^D^/2{F = 2) transition: 

/Yb = 688 358 979 309 307.6(1.4) Hz (1) 
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Figure 2: Results of absolute frequency measurements of the Yb+ optical frequency stan- 
dard at 688 THz, plotted as a function of measurement date (MJD: Modified Julian Date), 
covering the period from December 2000 to June 2006. The dotted straight line is the result 
of a weighted linear regression through the data points. The gray lines indicate the ±lcr 
range for the slope as determined from the regression. The intercepts of the gray lines are 
chosen such that they cross in the point determined by the weighted averages of frequency 
and measurement date. 



with a total relative uncertainty of 2.0 • 10~^^. This frequency value refers to a measurement 
at room temperature (296 K). Extrapolating to zero temperature (as it was done for the 
cesium clock) the frequency would have to be corrected for the AC Stark shift from blackbody 
radiation, i.e. increased by 0.37(5) Hz. 

From a weighted linear regression of the sequence of frequency measurements (cf. Fig. 
2) we obtain a slope of (—0.54 ± 0.97) Hz/yr, corresponding to a value for the fractional 
temporal variation of the frequency ratio dh\[fYh/ fcs) / dt — (—0.78 ± 1.40) • 10~^^ yr~"^, 
consistent with zero. This result can now be used to obtain new stringent limits on variations 
of a and /i. 

3. Limit on variations of the fine structure constant 

For the analysis we use a simple parametrization that includes only a minimum of as- 
sumptions. The electronic transition frequency is expressed as [4j [26] 



/ = const • Ry ■ F{a) (2) 

where Ry = niee'^ / (8eoh^) ~ 3.2898 ■ 10^^ Hz is the Rydberg frequency, appearing as the 
common atomic scaling factor. F(a) is a dimensionless function of a that takes relativistic 
contributions to the level energies into account [571 [H] • The numerical constant in front 
depends only on integer or half-integer quantum numbers characterizing the atomic structure 
and is independent of time. The relative temporal derivative of the frequency / can be 
written as: 

din/ dlnRy d\na d\nF 

— ; — = ; ^A- — ; — with A = — . (3) 

dt dt dt d\na ^ ' 

The first term d In Ry / dt represents a variation that would be common to all measurements of 
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electronic transition frequencies: a change of the numerical value of the Rydberg frequency^ 
The second term in Eq. 3 is specific to the atomic transition under study. The sensitivity 
factor A for small changes of a has been calculated for the relevant transitions by Dzuba, 
Flambaum et al. [281129]. 

To obtain a limit on dh\a/dt we combine the data from ^^^Yb"*" with results on the 
constancy of a transition frequency in ^^^Hg+, measured at NIST [3TJ[2S]. The transition 
5c?^'^6s — > 5rf^6s^ ^^5/2 of the mercury ion at 282 nm (1065 THz) with a natural 
linewidth of 1.9 Hz is studied. A sequence of measurements over the period from August 
2000 to November 2002 has resulted in a constraint on the fractional variation of /Hg//cs 
at the level of 7 • 10"^^ yr'^ gT]. From more recent data[25 with lower uncertainty a new 
limit dln(/Hg//cs)/di = (0.27 ± 0.78) ■ lO-^^ yr"! can be derived. 

The sensitivities of the ytterbium and mercury transition frequencies to changes of a are 
quite different [28l[29]: ^vb = 0.88 and Ang = —3.19. The change of sign of A between the 
two transitions reflects the fact that in Yb+ a 6s-electron is excited to the empty 5d-shell, 
while in Hg+ a hole is created in the filled 5d-shell when the electron is excited to 6s. The 
two measured drift rates together with Eq. 3 can now be used to calculate: 

^ = (-0.26 ± 0.39) • 10-15 yr'\ (4) 
(-0.55 ± 1.11) • 10"^^ yr-\ (5) 



dt 

We conclude that to within an uncertainty of 4 parts in 10^^ per year, the fine structure 
constant is constant in the present epoch. This result can be qualified as model-independent 
because it is only the simple expression Eq. 2 and the ab initio atomic structure calculations 
of the A-values that have been used in the interpretation of the experimental data. The 
constancy of Ry implies that there is no temporal drift between atomic gross structure 
frequencies and the cesium hyperfine splitting. This is of importance for metrology because 
it means that a cesium clock and an optical atomic clock can be regarded as fundamentally 
equivalent. 

Figure 3 shows how the two frequency drift rate measurements contribute to the con- 
straints on changes in a and in Ry. Shown is also the standard uncertainty ellipse within 
the central parallelogram. In comparison to the earlier analysis that we published in 2004 
|30| . based on the data available at the time, the Icr uncertainty in dlna/dt has decreased 
by a factor of five. Though over a very different time scale, the laboratory experiments have 
now reached the same sensitivity for dlna/dt as the analysis of quasar absorption spectra, 
if a linear time evolution of a is assumed. 



4. Limit on variations of the proton-to-electron mass ratio 

The ratio fi = mp/me is the second important dimensionless number whose possible 
variability was studied in astrophysical observations, especially in H2 and OH molecular 
spectra. In the laboratory, molecular spectroscopy could also provide a direct access to 
a variability of fi [42^, but has not yet reached the precision that has been achieved for 
electronic and hyperfine transitions in atoms. Data from the latter can, however, be used 
to derive a model-dependent limit on dlufi/dt along the lines described in Rcfs. [43,44]. 

^This statement refers to a dimensional quantity (frequency) in an SI unit (hertz). Since the SI second is 
fixed to a property of the '^^^Cs nucleus, possible variations of the quantity Ry are not purely electromagnetic 
in origin — as one would expect for an abstract physical quantity "Rydberg constant" - but are also related 
to the strong interaction. 
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Figure 3: Limits on temporal changes of the fine structure constant a and the Rydberg fre- 
quency Ry obtained from the measured drift rates of frequencies in Yb"*" and Hg+. Stripes 
mark the la uncertainty of the individual measurements and the resulting combined stan- 
dard uncertainty ellipse is shown within the central parallelogram. 



In analogy to equation 3 the cesium atomic hyperfine transition frequency may be 
parametrized as 

fc,=const- Rya^ -^-F' (a) (6) 
MS 

where F'(a) describes relativistic corrections and fJ-cs/fJ-B is the cesium nuclear magnetic 
moment in units of the Bohr magneton ^s- Consequently, the logarithmic change of the 
ratio of an optical electronic frequency /opt and fcs is given by 

c^ln(/opt//cs) dln(Mcs/Ms) , (Una 
dt = dt + " " ^> 

where A'^^ = dhiF^g/dlna ~ 0.83 [28 . Linear regression of Eq. 7 with the data from the 
optical frequencies in Hg+ and Yb+ gives the result 

dln(MCs/Ms) ,-,o_L01^ in-15..„-l 



dt 



1.3±2.1) • 10"'^ yr"\ (8) 



A very precise result has been published on the constancy of the ratio of ground state 
hyperfine frequencies in ^'^Rb with respect to ^^^Cs [20]: d\n{fcs/ fRh)/dt = (0.05 ± 0.53) ■ 
10~^^ yr~i. Since this frequency ratio scales like (/iCs/MRb)^"'*^ one may combine this 
experimental result with the limit on the variation of a (Eq. 4) in order to establish a limit 
on the temporal drift of the ratio of magnetic moments: 

rfln(MCs/MRb) _ ,n io_Ln -,n-i5 ,,,.-1 



dt 



= (0.18 ±0.56) • 10"^^ yr~^ (9) 



We can now derive model-dependent results on the proton (/-factor gp and on the proton- 
to-electron mass ratio by using the Schmidt model [45] that relates the nuclear magnetic 
moment to the magnetic moment of a single unpaired nucleon - a proton in the cases of ^^Rb 
and ^■^^Cs. The predictions of the Schmidt model agree with the actual values to within 
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about 25% for ^^Rb and to within 50% for -'^^■^Cs. An uncertainty due to the apphcation of 
this model is not included in the uncertainties of the limits presented in the following. From 
the Schmidt model expression /ics/MRb oc (10 — gp)/{gp + 2), the limit on the drift of the 
magnetic moments can be applied to constrain a possible drift of gp: 

= (-0.09 ±0.28)- 10-15 yr-i. (10) 

Using the Schmidt model expression for the cesium magnetic moment /ics = 7(10 — 
ffp)MJv/18 (with the nuclear magneton fi^) and the fact that ^iB/fJ'N — rnp/nie = ^ one can 
write the ratio of an optical electronic transition frequency to the cesium hyperfine frequency 
as 

/opt Fppt n , , 

fcs "^a^F^^ 10 -.g/ ^''^ 

Combining the limit on the drift rate of gp (Eq. 10) with the results on the constancy of 
the optical frequencies in Hg+ and Yb"*", a limit on the present temporal variation of the 
proton-to-electron mass ratio is obtained: 

^ = (-1.2±2.2).10-i'5yr-\ (12) 
at 

Compared to the previous evaluation |43j these limits to variations of the magnetic 
moments and mass ratio are now more strict by about a factor of three. Like in the case 
of a, the laboratory experiments reach the same sensitivity for dhifi/dt as the analysis of 
quasar absorption spectra under the assumption of a linear time evolution. We would like 
to point out, however, that the presented result on da/dt should be valued higher than the 
one on dfi/dt because the latter is model-dependent. A more detailed parametrization of 
nuclear magnetic moments in terms of quark masses has been published ,46J and can be 
used if more precision data from different nuclei becomes available. 

5. Conclusion 

The recent progress in the accuracy and reproducibility of optical frequency standards 
allows one to perform a sensitive laboratory search for temporal variations of fundamen- 
tal constants. So far all these experimental results are consistent with the constancy of 
constants. We have presented an update on limits to variations based on recent optical 
frequency measurements in single trapped ions. Coincidentally, the laboratory experiments 
over time scales At « 5 yr presently reach about the same sensitivity for changes of the fine 
structure constant Aa/At as the astrophysical observations in quasar absorption spectra 
over a time scale At « 10^" yr. The laboratory search will soon become more sensitive, as 
the precision of frequency standards continues to improve. In addition, a larger variety of 
systems is now being investigated so that it will be possible to perform tests of the consis- 
tency if a first observation of a variation in the laboratory is reported. The introduction of 
cooling techniques for small molecules will lead to more precise spectroscopy and to more 
sensitive tests for a variability of /i. In nuclear physics, a very promising system has recently 
been identified in the closely spaced lowest states of the ^^^Th nucleus [47j. This system has 
the potential to extend the scale for a search for changes in the strong and electromagnetic 
interactions by five orders of magnitude [48j . 
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